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extern "C" void 

sRGBColor( Color *pOut, Color *pln ) 
{ 

for ( int i = 0; 1 < 4; i++) 
{ 

float x = (*pln)[i]; 
If ( X < 0.03928f ) 

x = x/12.92f; 
else 

X = powf(((x + 0.055f)/1.055f), 2.4f); 
rpOut)[i] = x; 
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FIG. 3 
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extern "C" void 

LinearTosRGBColor( Color *pOut, Color *pln ) 
{ 

// coarse approximation: weigiited aritlimetic mean between 
// x^O.5 and x^O.375 approximates x'^(1/2.4) 
for ( int i = 0; i < 4; i++) 
{ 

float x = {*pln)[i]; 

float sqrtx = sqrtf(x); 

float sqrtSx = sqrtf(sqrtf(sqrtx)); 

float pow124 = 0.38f*sqrtx+0.62f*sqrtx/sqrt3x; 

if ( x < 0.00304f ) 

x = 12.92f*x; 

else 

x = 1.055f*pow124-0.055f; 
rpOut)[i] = x; 

} 

} 
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extern "C" void 

LinearTosRGBColor( Color *pOut, Color *pln ) 
{ 

// finer approximation that avoids taking 3 successive 
// square roots: apply one round of N-R to guess cube 
//rootofx*sqrt(sqrt(x)) 
for ( int i = 0; i < 4; i++) 

{ 

float X = (*pln)[i]; 



float sqrtx = sqrtf(x); 

float sqrt2x = sqrtf(sqrtx); 

float appx = 0.78f*sqrtx+0.22rsqrt2x; 

float num = x*sqrt2x; 

float cuberoot = (2*appx+(x*sqrt2x)/(appx*appx))/3.0f - 0.00025f; 
if ( x < 0.00304f ) 

x = 12.92f*x; 

else 

X = 1.055rcuberoot-0.055f; 
(*pOut)[i] = x; 

} 
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const declspec(align(16)) m128 Const039 = _mm_set1_ps( 

0.03928f ); 

const ^declspec(align(16)) _m128 Constlnv1292 =_nnm_set1_ps( 

1.0f/12.92f ); 

const _declspec(align(1 6)) _m128 Const055 = _mm_set1_ps( 0.055f 
); 

const ^declspec(align(16)) m128 Constlnv1055 = _mm_set1_ps( 

1.0f/1.055f ); 

const _declspec(align(16)) _m128 Const1285 = _mm_set1_ps( 
1.285f); 

const _declspec(align(16)) _m128 Const0285 = _mm_set1_ps{ 
0.285f ); 
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extern "C" void 

sRGBColor( Color *pOut, Color *pln ) 
{ 

// SIMD: compute BOTH answers and compose output using mask 
_m128 ansBelowDelta =_mm_mul_ps( *(_m128 *) pin, Constlnv1292 ); 

m128 X = _mm_mul_ps( Constlnv1055, _mm_add_ps( *( m128 *) pin, 

Const055 ) ); 

m128 sqrx = _mm_mul_ps( x, x ); 

ml 28 invsqrx = _mm_rcp_ps( sqrx ); 

m128 invsqrtx = _mm_rsqrt_ps( x ); 

ml 28 ansAboveDelta = _mm_div_ps( Const1285, 

_mm_mul_ps( invsqrx, _mm_add_ps( Const0285, invsqrtx ) ) ); 
_m128 TruelfLTDelta = _mm_cmplt_ps( *(_m128 *) pin, Const039 ); 
*{ ml 28 *) pOut = _mm_or_ps( _mm_and_ps( TruelfLTDelta, ansBelowDelta 

), 

_mm_andnot_ps( TruelfLTDelta, ansAboveDelta ) ); 
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const_decispec(align(16)) _m128 CONST00304 = __mm_set1_ps( 0.00304f ); 
const _declspec(aiign(16)) _m128 CONST1292 = _mm_set1 js( 12.92f ); 

const declspec(align(16)) m128 CONST055 = _mm_set1_ps( 0.055f ); 

const _dec!spec(align(16)) _m128 CONST1055 = _nnm_set1_ps( 1 .055f ); 

const declspec(aiign{16)) m128 CONST078 =_mm_set1_ps( 0.78f ); 

const _declspec(align{16)) _m128 C0NST1 m078 = _mm_set1_ps( 1 .Of-078f ); 

const _declspec(align(16)) _m128 CONST38 = _mm_set1_ps( 0.38f ); 
const declspec{a!ign(16)) m128 CONST1m38 = _mm_set1_ps( 1 .0f-0.38f ); 

extern "C" void 

LinearTosRGBColor( Color *pOut, Color *pln ) 
{ 

m128 ansBelowDelta = _nnm_nnul_ps( *( m128 *) pin, CONST1292 ); 

nn128 sqrtx = _mm_sqrt_ps( *( m128 *) pin ); 

ml 28 sqrt3x = _mm_sqrt_ps( _mnn_sqrt_ps( sqrtx ) ); 

_nn128 pow124 = _mm__add_ps( __mm_muLps( CONST38, sqrtx ), 

_mm_div_ps( _mm_muLps( CONST1m38, sqrtx ), sqrtSx ) ); 

nn128 ansAboveDelta = _nnnn_sub_ps( _nnm_muLps( CONST1055, pow124 ), CONST055 ); 

m128 TruelfLTDelta = _mm_cmplt_ps( *( m128 *) pin, CONST00304 ); 

*(_m128 *) pOut = __mm_or_ps(_mm„and_ps( TruelfLTDelta, ansBelowDelta ), 
_mm_andnot_ps( TruelfLTDelta, ansAboveDelta ) ); 

} 
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const declspec(align(16)) m128 Magic00304 = _mnn__set1_ps( 0.00304f ); 

const _declspec(align(16)) _m128 Magic1292 = _mm_set1 jds( 12.92f ); 
const_cleclspec(align(16))__m128 Magic055 =_nrim_set1__ps( 0.055f ); 
const_declspec(align(16))_m128 Magic1055 =_mm_set1_ps( 1.055f ); 

const declspec{align(16)) m128 MagiclnvS = _mm_set1_ps( 1 .0f/3.0f ); 

const declspec(align(16)) m128 MagicFudge = _!nm__set1 jds( 0.00025f ); 

const _declspec(align(16)) _m128 Magic078 = _mm_set1_ps( 0 J8f ); 
const_declspec(aiign(16))_m128 Magic1nn078 = _mm_set1_ps( 1.0f-0.78f ); 

const declspec(aiign(16)) m128 Magic38 = __mm_set1_ps( 0.38f ); 

const declspec(align(16)) m128 Magic1m38 = _min_set1_ps( 1.0f-0.38f ); 

extern "C" void 

LinearTosRGBColor{ Color *pOut, Color *pln ) 
{ 

m128 ansBeiowDelta = _mm_mul_ps( *( m128 *) pin, Magic1292 ); 

m128 sqrtx = _mm_sqrt_ps( *( m128 *) pin ); 

m128 sqrt2x = _mm_sqrtjDs( sqrtx ); 

m128 appx = _mm_add_ps( _nnm_mul_ps( Maglc078, sqrtx ), 

_mm_nnul_ps( Magic! m078, sqrt2x ) ); 

m128 cuberoot = _mm_sub__ps( 

_mm_mul_ps( Magiclnv3, 
_mm_add_ps( _nnnn_add_ps( appx, appx ), 

_nnm_div_ps( _mnn_mul_ps( *( m128 *) pin, sqrt2x ), 

_mm_mul_ps( appx, appx ) 

) 

) 

), MagicFudge 

): 

_m128 ansAboveDelta = _mm_sub_ps( _mm_mul j)s( Magic1055, cuberoot ), Magic055 ); 
_m128 TruelfLTDeita = _mm_cmplt_ps( *(_m128 *) pin, Magic00304 ); 

*( nn128 *) pOut = _mm_or_ps( _nnnn_and_ps( TruelfLTDelta, ansBelowDelta ), 

_nnm_andnot_ps( TruelfLTDelta, ansAboveDelta ) ); 

} 
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